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I. INTRODUCTION 


This research aimed to improve the accuracy of critical 
wave steepness and breaking wave steepness estimation. 
The breaking wave steepness equation is derived from the 
wave amplitude function obtained through the integration 
of KFSBC over time (Hutahaean, 2024). Critical wave 
steepness is determined by identifying the maximum wave 
amplitude in the dispersion equation, where the dispersion 
equation results from applying the wave amplitude function 
to the Euler conservation momentum equation. Therefore, 
both critical wave steepness and breaking wave steepness 
are derived from the wave amplitude function. 


Hutahaean (2024) conducted integration with zeroth-order 
accuracy, yielding a wave amplitude function that produced 
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Abstract— The relationship between wave amplitude and 
other wave constants is derived through the integration of the 
Kinematic Free Surface Boundary Condition (KFSBC) over 
time, yielding the wave amplitude function. This function 
presents the characteristics of wave breaking to determine the 
breaker length index—the ratio of breaking wave height to 
breaking wave length. Previous research achieved integration 
at a zero-order accuracy level, assuming small wave 
amplitude and long wave, resulting in a notably large critical 
breaker length index and wave steepness. In this research, 
KFSBC integration was advanced to second and third-order 
that 
significantly reduces the critical breaker length index and 


accuracies, yielding a wave amplitude function 


compared to zero-order integration. 


the 


a comprehensive wave 


wave steepness 


Subsequently, utilizing third-order wave amplitude 


function, transformation model 
incorporating shoaling, breaking, refraction, and diffraction 


was developed, demonstrating robust model performance. 


significant critical and breaking wave steepness. This 
research extended the integration to second and third-order 
accuracy levels, examining resulting critical wave steepness 
and breaking wave steepness. 


Additionally, Hutahaean (2024) developed a water wave 
transformation model based on the wave amplitude 
function. This research further refined the model with a new 
wave amplitude function, investigating breaking parameters 
such as breaking wave height and depth. The findings 
demonstrate that higher integration accuracy improves the 
wave amplitude function's effectiveness in predicting these 
parameters, validating the chosen integration method's 
efficacy. 
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Il. COMPLETE SOLUTION OF THE 
LAPLACE EQUATION 


The complete solution of the weighted Laplace Equation 
using the separation variable method, Hutahaean (2024) in 
deep water where the bottom slope has no effect, is 


(x,z,t) = 

G (cosk,x + sink,x) coshk,(h+z)sinot .....(1) 
(x, z,t) : velocity potential 

G : wave constant 


k, : wave number on the horizontal-x axis 


k, : wave number on the vertical-z 


po 
“We 


k is the general wave number 


Yx and y; is weighting coefficient in weighted Taylor series, 
see section (3). 


20. . : 
= is angular frequency, where T is wave period. 


h : water depth 


Equation (1) is carried out at the characteristic point where 
cosk,x = sink,x, 


o(x,z,t) = 2G cosk,xcoshk,(h+z)sinot .(2) 


IW. WEIGHTED KINEMATIC FREE SURFACE 

BOUNDARY CONDITION 
a. Weighted Taylor series 
The Weighted Taylor series refers to a modified Taylor 
series truncated to first-order terms, where each first-order 
term incorporates a specific weighting coefficient in place 
of higher-order terms (Hutahaean, 2023). 
Weighted Taylor series for function f = f (x, t), 

df 


df 
f(x + 6x,t + 6t) = f(x,t) + Ye2ôt r + YxOX T 


Weighted Taylor series for function f = f (x, z, t), 
f(x + 6x,z + 6z,t + ôt) 
df 
= fz, t) + ¥r30¢ mr + 
a a 

yrOx b+ yz E ...... (4) 
Yt3> Yt2 Yx and y, are weighting coefficients. The basic 
values of the weighting coefficient are, y,3 = 3.0 , Yt2 = 
2.0 , y,=1.0 and y,=1. The corrected weighting 


coefficient values as a function of the optimization 
coefficient £ are presented in Table (1) as follows. 
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Table (1). Corrected weighting coefficients values. 


E Yt2 Yt,3 Yx Yz 
0.010 1.9998 | 3.00465 | 0.99879 | 1.01093 
0.011 | 1.99975 | 3.00563 | 0.99854 | 1.01325 
0.012 | 1.99971 | 3.00671 | 0.99826 | 1.0158 
0.013 | 1.99966 | 3.00788 | 0.99795 | 1.01858 
0.014 | 1.99960 | 3.00915 | 0.99763 | 1.02159 
0.015 | 1.99954 | 3.01052 | 0.99727 | 1.02484 
0.016 | 1.99948 | 3.01198 0.9969 | 1.02832 
0.017 | 1.99941 | 3.01355 | 0.99649 | 1.03205 
0.018 | 1.99934 | 3.01521 | 0.99607 | 1.03601 
0.019 | 1.99926 | 3.01697 | 0.99561 | 1.04022 
0.020 | 1.99918 | 3.01883 | 0.99514 | 1.04468 

b. Weighted Kinematic Free Surface Boundary 


Condition (KFSBC). 

For the water surface elevation equation n = n(x,t), the 
weighting Taylor series is, 

an 

dx 

The first term on the right side is moved to the left, and the 
Equation is divided by ôt, 


d 
n(x + ôx, t + ôt) = n(x,t) + yeaðt E + yxôx 


n(x + dx,t + ôt) — n(x, t) dn ôx dn 

ôt ~The Tt 
As ôt approaches zero, the left side of the equation 
represents the speed of change in water surface elevation, 
which corresponds to the total vertical velocity of surface 
water particles. 


Wy = Yt,2 = + Yxun a ...(5) 

This equation is denoted as KFSBC. Here, w, represents the 
vertical velocity of surface water particles, and uy denotes 
the horizontal velocity at the water surface. 


Equation (5) can thus be expressed as the equation 
describing changes in water level elevation, namely, 


dn dn 
= = Wy — Yup — (6 
Yt,2 dt n Vx n ax (6) 


IV. KFSBC INTEGRATION WITH 0TH ORDER 
ACCURACY. 

Hutahaean (2024) utilized zeroth-order integration to 

establish the relationship between wave amplitude and 

wave constants G, such as wave number k and angular 

frequenc o, termed as the wave amplitude function. The 

subsequent section will reiterate the steps involved in the 
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zeroth-order integration process to outline the formulation 
of the equation for the wave amplitude function. 


1.1. The formulation of wave amplitude function. 
Using equation (2), the vertical water particle velocity is: 
db 


w= ca S —2 Gk, cosk,x sinh k;(h + z) sin ot 


Surface vertical water particle velocity is 
W, = —2 Gk, cosk,x sinh k,(h + n(x, t)) sinot 


Meanwhile, the horizontal water particle velocity is 
d 
u= -— Z? = 2Gk,, sink,x coshk,(h + z) sin øt 


Surface horizontal water particle velocity is 
Un = 2Gk, sink,x cosh k,(h +n(x, t)) sin ot 


The substitution of (7) and (8) number (6), 
dn 


Vezme. dt 
—2G k, cos k, sinh k; (h + n(x, t)) sinot 


_ dn 
—2y,Gk,, sin ky cosh kz (h + n(x, t)) sin ot ax 


To derive the equation for water surface elevation, integrate 
this equation over time t, 


Ye2n(x,t) = 

—2Gk, cos kx | sinh k,(h +n, t)) sin ot dt 
—2yyGk, sin k,x f coshk,(h + 

n(x, t)) sinot = dt ..(9) 


The water surface elevation equation obtained depends on 
the integration solution's accuracy level. Theoretical 
accuracy ranges from order 0 to order n, determined by the 


a"n. : R : 
presence of —— in the integration solution. 
at” 


; : aon . i 
Order 0 Da assumes an ie in the results, while order 


0 


1 implies = , etc. This research was 


= order 2 involves “2 aa 


limited to wen 3 accuracy. 
In zero-order integration, the assumption is made that in 


deep water and with small wave amplitudes |sinh kh (1 + 


1a) and |cosh k,h (1 + 10) a7] fluctuates very little 


with time, allowing it to be treated as constant and excluded 
from integration. 


Ye2n(x,t) = 
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—2Gk, sinh k,(h +n(x, t)) cos k,.x | sin ot dt 
—2y,Gk, cosh k; (h 
. dn f 
+n(x, t)) sin k,x = | sin ot dt 


This integration can be completed by integrating the 
sinusoidal function alone. After integrating, the 


characteristic point property is worked out, namely 
$ . . k 
cos kx =sink,x and substituted with k, =-= and 


VYx 


k 
V2 


2Gk 
Ye2n(%,t) = a coshk z(h +(x, t)) 


(= kz it n(x, OD), + hex =) cosk,x cos ot 


kz = 


In the integrated results, the absence of the a term indicates 
a zero-order integration. 

In deep water, where tanh k,(h + n(x, t)) = 1, k(h + 
n(x,t)) = On. @ is referred to as deep water coefficient 
with 0 > 1, for instance 0 = 1.75. 


cosh 61 


2Gk 
n(x, t) = 


on On 
(Z Yx cos kx cos ot 
VYz Weg az) 


As a periodic function, hence 


2Gk tanh 67 dn 
A = — hé ( Af =) 
ay cosh Om ie + V/V Ax 


t,2 


A is wave amplitude. 


Water surface elevation equation is 


n(x,t) = Acosk,x cos ot ...(10) 
At the characteristic point of space and time, 

dan — kA 

ax 2 


The wave amplitude function is 


2Gk tanh 67 kA 
loon cosh On (FES — 5 . (11) 


This equation is the wave amplitude function equation 
resulting from zero order integration. 


1.2. Oth order breaking characteristics. 
In (11), breaking occurs when 
tanhðm kA | i 
Vz 2 
Hp = 2 tanh 0r 
ip ORE (12) 
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Equation (12) applies universally across all wave periods, 
accommodating the specific wave height relevant to each 
period. 


1.3. 0" order dispersion equation 

The dispersion equation is formulated based on surface- 
weighted Euler momentum conservation, excluding the 
convective acceleration term. 


dun _ an 


Yt Ge TTI ax 


With u, from (8) and 7 and (10), the following is obtained 


du, ; 
Yt3 Er 2Yt3Gkx0 sin kx cosh k,(h 


+n(x, t)) cos ot 
dn ; 
—g ax gk,Asink,x cos ot 
Substitute these two equations into (13), 
2¥¢3Go0 cosh Om = gA 
Wave amplitude on the right side is substituted with (11), 


TE né 2gGk ne = On A 
Yt3G0 cosh Or = cosh 0r — — 
G OYt2 VYz 2 

A tanh 8 
= k? = = + Vt 2Yt30° = 0 PEA (13) 


This equation is the dispersion equation resulting from zero 
order integration. 


The steps in formulating the wave amplitude function are: 

a. Completing the integration of KFSBC 

b.The integration results are collected to form an 
Equation:n(x,t) = (..... ) cos k,x cos ot 

c. Defining wave amplitude function A =(...... ), obtaining 
n(x,t) = Acosk,xcosot 


. . ; a E 
From this Equation, we obtain the Equation of oe which is 


then substituted into the wave amplitude function equation, 
to obtain the final Equation of the wave amplitude function. 


V.  KFSBC INTEGRATION WITH 2ND ORDER 
ACCURACY. 
The integration will be completed until there is an element 


an i : : : : : : 
az" Integration is carried out using the integral inversion 
method (Hutahaean (2010). 


a. f sinhk,(h + n(x, t)) sinot dt 


Is defined as, 


f(t) = sinh k,(h + n(x, t)) cos ot 
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df 
FT al sinhk,(h + n(x,t)) sinot + 


d 
k, coshk, (h +n(x, t)) cos ot z 


This differential equation is integrated over time t, 
[ a = -0 Í sinh k,(h + n(x,t)) sinot dt 


dn 
+k, Í cosh k,(h +n(x, t)) cos ota dt 


Integration of the left side produces f(t), substitution 
f (t)and the rearranged Equation is obtained 


| sinh k,(h + n(x, t)) sinot dt = 
1 
— 5 sinh k,(h + n(x, t)) cos ot 


k d 
+ = | cosh k,(h + n(x, t)) cos ot dt 


In second-order integration, the fluctuations 
t)\ a . . 
| coshk,h (1 + 10) 2| with respect to time are 


significant enough that the integration of the two terms on 
the right side necessitates the use of the integral inversion 
method. 


dn 
| cosh k,(h + n(x, t)) cos ot dt =? 
Is defined as, 


_ dn 
f(t) = cosh k,(h +(x, t)) sin otar 


2 


df dn 
u~ k, sinh k,(h +n(x, t)) sin ot (=) 


dn 
+o cosh k, (h + n(x, t)) cosot at 
2 


„an 
+ cosh k,(h +n(x, t)) sinot ap 


Is integrated into, 


fa- 
. . any? 
i k, sinh k,(h +n(x, t)) sin ot (=) dt 


dn 
+a f cosh k,(h +n(x, t)) cos ota dt 


d*n 
+ | cosh k,(h +n(x, t)) sin otra dt 


The left side is substituted with f (t) and rearranged 
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dn 
| cosh k,(h +n(x, t)) cos ota dt = 


1 _ diy 
5 cosh k,(h + n(x, t)) sin ot ar 


ke [ l m 
-= | sinh k,(h + n(x, t)) sin ot (=) a 


1 Ge 
— =Í cosh k,(h +n(x, t)) sin ota dt 
In this 2nd order integration, the assumption is made that 


fluctuations over time t from |sinhk,h (1 + 


Foley 
h dt 


small and can be considered constant that it can be excluded 


is very 


and |cosh kz n(4 + aD) 2a FE 


from the integral. 


dn 
| cosh k,(h +n(x, t)) cos ota dt = 


1 _ dn 
5 cosh k,(h + n(x, t)) sin otar 


kz an? f 
-sinh k,(h + n(x, t)) (=) | sinot dt 


1 d2n 
= 5 cosh k,(h + n(x, Dl sinot dt 


The final integration is completed by simply integrating the 
sinusoidal function, obtaining the followings. 


dn 
| cosh k,(h +n(x, t)) cos ota dt = 


1 _ diy 
5 cosh k,(h +n(x, t)) sinot ar 


2 


k d 
+ <= sinh k,(h + n(x, t)) cosot (=) 


1 €?n 
+ zz cosh k z(h +n(x, t)) cos ot —> a? 
Or, 


k dn 
=| cosh k,(h + n(x, t)) cos ote dt = 


kz _ dn 
G2 cosh k,(h + n(x, t)) sinot ar 


2 


k2 d 
+ <3 sinh k, (h +n(x, t)) cos ot (=) 


k, d?n 
+ cosh k z(h + n(x, t)) cos ot > 3 


Since k, (h + n(x, t)) = r and calculating the integration 


outcomes at the specific time point, cos ot = sinot 
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k dn 

=| cosh k,(h + n(x, t)) cos ota dt = 
1 k,dn_ k? d 
— cosh 61 aa + — tanh 0r (=) 
o odt o? dt 


k, d2n 
++; a ar cos ot 


2 


The results of this integration are substituted into (14), 


cosh 61 


f sinn k,(h + n(x, t)) sinot dt = — = 


k,dn_ k2 an’? 
tanh 0r — | —— + — zta nh on (= =) 
odt o? dt 


k, @°n 
o* dt? 


cosot........(15) 


i : a? : 
Based on the integration results, there are terms a that is 


; dn\? ae ck : 
equivalent to (=) , classifying this integration as second- 


order accurate. 


b. fcoshk,(h + (x, t)) sin ot 2 dt 


Is defined as, 


dn 

f(t) = cosh k,(h + n(x, t)) cos ots 
_ dn 
= —ocosh k,(h + n(x, t)) sin ota 


dn d 
+k, sinh k z(h +n(x, t)) cos ot = = 


2 


+coshk,(h + (x, t)) a! 
coshk, n(x,t)) cosa TA 


Multiplied by dt is integrated into, 
„dn 
| cosh k,(h +n(x, t)) sin ote dt = 


1 dn 
— 5 cosh k, (h + n(x, t)) cos ot ax 


k dnd 
+f sinh k,(h +n(x, t)) cos a dt 


+f hk,(h + (x, t)) cos ot Ti 
COS NX, CcOSO dtdx 


In this second order integration, the assumption is made that 
: ; t)\ ana 

the fluctuation |sinh k,h (1 + a: 2) ga 

time t is quite small that it can be considered constant and 


| with respect to 


can be excluded from the integral. 
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with 


Likewise, the fluctuation I|cosh k,h (1 + ul xa?) RTE 


respect to time t is considered too small that it can be 


excluded from the integral. 
_ dy 
| cosh k,(h +n(x, t)) sin ote dt = 
1 d 
— 5 cosh k, (h +n(x, i cos ote 


kz 
+ sinh k z(h +n(x, n) cosot dt 


1 d?n 
+z cosh k,(h +n(x, ae ot dt 


This integration equation can be solved by integrating 
fcosot dt directly. Once integrated and substituting 
k,(h+n(x,t)) = On, the equation is evaluated at the 
characteristic time point where cos ot = sin ot, yielding: 


f coshk,(h + n(x, t)) sin ot dt = 


—+cosh 6 
o 
dn kz nan 
=— — — tanh 0x — — 
(2 anaon at ax 
i t...(16 
- z gray | OS (16) 
Regarding the presence of 2 — Z an nd = — z then the result of 


the integration is called i i ader level of accuracy. 


Equation-Equation (15) and (16) are substituted to (9), 


Gk, 
n(x, t) = cosh 61 
OY t,2 


( tanh 61 


k2A? 
+ (raa = I tanh on) cos ot cosk,x 


Gk d k dnd 
gp ee osh On (2 ——tanh jp 
dx o 


OY t,2 


cos ot cos kx 


As a periodic function, therefore 


ee Gk, 
OY t,2 


cosh 61 ( tanh 61 


k2A? 
+ (raa = i tanh on) 
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Vx Gk dn kz dn dn 
—— cosh 0 aoe nh 62 
j OYt,2 — (F at dx 
1 in 
o dtdx 
obtaining 


n(x,t) = A cosotcosk,.x 
Within the characteristic point of 


dn — kA 
dax 2 
dn oA 
dt 2 
d*n ok, A 
dtdx 2 
Is the differential substitution of y(x,t) and k,= 
-= and k, = = 
V¥x 7 Vya 
A= cosh 61 


OY t,2 


Cave 


e tanh 67 
k? A? 
t 4 


This equation is the wave amplitude function resulting 
from integration with 2nd order accuracy. 


5.2. Breaking 2nd order characteristics 


Breaking occurs when 


tanh 61 r (- 1) kA 
VYz z 
1 tanh 67 
( + 2) ——— k*A* =0 
VeV¥z V¥z) 4 


: 2 F ; H 
Is substituted to k = = and sinusoidal wave A = a 


Samer -aef - 
a a 
=0 . (18) 


Using equation (18), “known as the breaker length index ae 
b 


can be calculated. This equation does not include 


parameters for wave period or wave amplitude. Therefore, 
H : : : : : 
obtained is valid for all wave periods, provided the 
b 


appropriate wave amplitude is considered. 
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5.3. 2nd order Dispersion Equation. 

In the same way as the formulation of the dispersion 
equation in (4.3) where in this case the wave amplitude 
function (17) is used, the dispersion equation obtained is 


( 1 ig ve) tanh 67 


1 
gA*k3 — C- — 1) gAk? 
Z 


YzVYz VYz) 4 
tanh 67 
a a gk + Yt,2¥t,39 = 0 sa (19) 


The equation is a third-degree polynomial that can be solved 
using the Newton-Raphson iteration method. Initially, 
terms containing k? are disregarded to simplify the equation 
to a second-degree polynomial. Once k is found from this 
simplified equation, it is used in the complete set of 
equations to calculate k accurately. 


VI. KFSBC INTEGRATION WITH 3RD ORDER 
ACCURACY 


In this section the integration process is not discussed in its 
entirety, just an example. In second order integration there 
are integration results, 


dn 
| cosh k,(h +n(x, t)) cos ota dt = 


1 _ diy 
5 cosh k,(h + n(x, t)) sin ota 


ka [ i am 
== Í sinh k,(h + n(x, t)) sinot (=) dt 


1 dn 
= =Í cosh k,(h +(x, t)) sin Otz dt 


In second-order integration, the two integrations on the right 
side are solved by integrating only the sinusoidal elements. 
However, in third-order integration, the two right-hand side 
integrations are solved using the integral inversion method. 
For example, 


k f . any 
=| sinh k,(h + n(x, t)) sinot (=) dt 
2 


f(t) = sinh k,(h + n(x, t)) cos ot E» 


df any? 
u” k, cosh k,(h + n(x, t)) cos ot (=) 


d 
—o sinh k,(h +n(x, t)) sinot (2) 


an dn 
+2 sinh k,(h + n(x, t)) cos ota az 
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This equation is multiplied with dt and is integrated into, 


fu- 


dn 
kz | cosh k,(h + n(x, t)) cos ot (=) dt 


. i a 
= | sinh k,(h + n(x, t)) sin ot (=) a 


an a’n 
+2 f sinn k,(h + n(x, t)) cos ota az dt 


Substitute f(t) on the left side and is moved to the right, 
while the second term on the right side is moved to the left, 


any 
| sinh k,(h + n(x, t)) sinot (=) dt = 


i any? 
=~ sinh k,(h + n(x, t)) cosot (=) 


dt 
d?n 


iA | nein ene ace = — cae 
= | sin z NX, cos o at ate 


k dan? 
+72 Í cosh k,(h + n(x, t)) cos ot (=) dt 


The assumption is that fluctuations with time from: 


[cosh k,h (1 + 102) (2) and |sinh k,h (1 + 


dt 


(x,t)\ dn d2n| . : 
n : ay S is very small and can be considered constant, 


that it can be excluded from integration, and what is 


integrated into, is only the sinusoidal element. Obtaining, 
k any? 
=| sinh k,(h +n(x, t)) sinot (=) dt = 
o 


kz dn 
— z sinh k,(h +n(x, t)) cos ot (=) 


3 


k2 d 
+ zs cosh k,(h +n(x, t)) sin ot (=) 


dt 
kz an dn 
+ z3 sinh k,(h + n(x, t)) cos otar az 


After all the integrations are completed then add them up 
and carry out the process as in 2nd order integration, we get, 


2Gkcosh 0r 
n(x, t) = ———— a(k, A) cos kx cos ot 
(k, A) tanh 01 ( 3 5) kA 
a š = SS 
VYz 2y, 2 
k? A? 
— tanh 0r 
2 SY NYz 
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é T 
YF Vez 


Where the wave amplitude function is 


2Gkcosh On 


A= CEA ....(20) 
Breaking occurs when, 
tanh On 3 3 (kA)? 
“~—+(—-5) — ——— tanh r 
VYz 2y, 2 2) YxVYz 
"e ae =0 (21) 
oop p es 


From this equation, the kA value can be calculated. It is a 
third-degree polynomial equation, which can be solved 
using the Newton-Raphson iteration method. Initially, the 
kA value is estimated by disregarding the terms. After 
obtaining an approximate kA value using this simplified 
equation, it can then be refined using the complete set of 
equations. For example, 


kA=A 
Considering k = = and H = 2A, 
H, A 
Ll, T 


Dispersion equation with wave amplitude function obtained 
from 3rd order integration is, 


1 1\A3\, tanhOrA? , 71 3), 
( = ) kt k +( = ) Ak 
Y? Yz/ 8 2 SYN Yz 2y, 2 
tanh 0r i ¥t3Vt3 0° 


+ — a 
VYz g 


=0 ...(22) 


This equation is a 4th-degree polynomial, which can be 
solved using the Newton-Raphson iteration method. The 
initial iteration value is obtained by disregarding the terms 
with powers of 4 and 3. 


VII. WAVELENGTH, CRITICAL WAVE 

STEEPNESS AND BREAKER LENGTH INDEX 
In this section, a comparative research of wavelength, 
maximum wave height, and critical wave steepness is 
conducted. The research utilizes the dispersion equation 
with an optimization coefficient £ of 0.01 for the weighting 
coefficient. 


In the dispersion equation derived from zero-order 
integration using the wave amplitude function, the 
maximum wave height corresponds to the wave height at 
which the determinant value of equation (13) becomes zero. 
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For the dispersion equation (19) obtained from second- 
order integration using the wave amplitude function, using 
an excessively large wave height results in a negative wave 
number. Therefore, the maximum wave height in equation 
(19) is the largest wave height that does not produce a 
negative wave number. 


Similarly, in the dispersion equation (22) resulting from 
third-order integration with the wave amplitude function, 
using a wave height that is too large also yields a negative 
wave number. Thus, the maximum wave height in equation 
(22) is the highest wave height before encountering a 
negative wave number. 

There are two types of wave numbers: the horizontal wave 
number k, and vertical wave number k,, where 


k k 
kx = Ta and kz = IE 

Therefore, two wavelengths are, 

Horizontal wavelength : L, = = = = Vex = LY Vy 


Vertical wavelength : L; = = = = VYyz = LV Vz 


Table (2) presents the maximum wave height of (13) 


Table (2). Maximum wave height from (13) 


T Hmax Iy Iz Ba 

(sec) (m) (m) (m) Lx 
6 1453 | 4.64 | 4.699 | 0.313 
7 1.977 | 6401 | 6481 | 0.309 
8 2.583 | 8.272 | 8.377 | 0.312 
9 3.269 | 10.487 | 10.619 | 0.312 


10 4.036 12.914 | 13.077 0.313 


11 4.883 15.708 | 15.906 0.311 
12 5.812 18.56 18.795 0.313 
13 6.821 21.79 22.064 0.313 
14 7.911 25.188 | 25.506 0.314 
15 9.081 29.056 | 29.422 0.313 


Table (3). Maximum wave height from (19) 


T Hoos ig Lz Hmax 
(sec) (m) (m) (m) | Ly 
6 1.565 6.169 6.246 0.253 
7 2.13 8.41 8.516 0.253 
8 2.782 10.988 11.126 0.253 
9 3.522 13.76 13.934 0.256 
10 4.348 17.027 17.241 0.255 
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11 5.261 20.618 | 20.878 0.255 
12 6.261 24.542 | 24.852 0.255 
13 7.348 28.799 | 29.162 0.255 
14 8.522 33.389 33.81 0.255 
15 9.783 38.31 38.793 0.255 


Based on the data presented in Tables 2, 3, and 4, it is 

evident that the difference between the horizontal 

wavelength L, and the vertical wavelength L,is minimal. 
Table (4). Wave height maximum from (22) 


T | Bie h a B 
(ec) | m | (m | mM | LE 
6 1.53 6.201 6.279 0.247 
7 2.082 8.468 8.575 0.246 
8 2.42 11.023 11.162 0.247 
9 3.444 13.839 14.013 0.249 


10 4.252 17.07 17.285 0.249 


11 5.146 20.474 | 20.733 0.251 
12 6.124 24.437 | 24.745 0.251 
13 7.186 28.835 | 29.199 0.249 
14 8.336 33.015 | 33.432 0.252 
15 9.568 38.292 | 38.775 0.25 


As the integration order increases, the critical wave 
steepness decreases. Specifically, the difference between 
the critical wave steepness values derived from equations 
(19) and (22) is relatively small, typically in the third 
decimal place. This suggests that increasing the integration 
order to 4th order would yield a critical wave steepness not 
significantly different from that obtained at 3rd order. In 
general, the critical wave steepness is estimated to be 0.250. 


Toffoli et al. (2010) propose a critical wave steepness of 
0.170, while suggesting it could reach up to 0.20. It is noted 
that the wave steepness calculated at 3rd order accuracy 
closely approximates the criteria proposed by Toffoli et al. 


‘ H 
The comparison of the breaker wave steepness 7 from 
b 


equations (12), (18), and (21), using parameters 0 = 1.75 
and €= 0.01, progresses from the lowest to highest 
integration order. 


Dari (12), ae = 0.629 
b 

Dari (18), <b = 0.446 
b 


Dari (21), ae = 0.437 
b 
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The critical wave steepness in deep water differs from the 


: . H 4 
wave steepness at the breaking point, =, where breaking 
b 


waves concentrate a significant amount of energy. Among 


H ) ; F 
the three values, the smallest is obtained from equation 
b 


(21), which corresponds to 3rd order integration. The 
difference between “tfrom equation (18) and “from 
b b 
equation (21) is negligible, indicating that increasing the 
integration accuracy will not lead to substantial changes. 


: ‘ Hp. 
In conclusion, the breaking wave steepness zis generally 
b 


estimated to be 0.437 


VII. WAVE TRANSFORMATION MODEL 
The wave amplitude function equation, which establishes 
the relationship among various wave constants, serves as a 
foundational equation for developing wave transformation 
models (Hutahaean, 2024). By deriving a new wave 
amplitude function, different wave transformation 
phenomena can be modeled. 
A 


y 


č 
a) N 
x 


Fig (1)Axis system for wave transformation modeling. 


The wave transformation model is constructed within the 
coordinate system depicted in Figure (1), where the wave 
propagates along the & axis in the (x,y) plane, forming an 
angle a with the horizontal x. 

8.1 Shoaling-breaking model. 

Shoaling breaking equations were developed for waves 
moving on the €é-axis. The wave amplitude function 
resulting from 3rd order integration, namely (20), is written 


as, 
A= Gk A(k, A) w1.(23) 
Ak, A) = 2cosh 01 (= On A ( 3 5) k ) 
i OYt,2 VYz 2y; 2 
2cosh ( k? 4? ne 
— anh 0r 
OY t,2 2) VeVVe 
A ( 1 1 ) =) 
Y Yz/ 8 
Equation (23) is differentiable to axis-€, 
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- =(65 +k E) Atk. A) (24) 


Where en the law of conservation of wave number 


; i da 
which will be shown in the next section, then aE = 


The conservation of energy equation (Hutahaean (2024) 
applies, 
dk 


This equation is written as, 
dG Gdk 
ae 2aE 
Substituted to (24), obtaining 
dA _ Gdk 


a Lue AKAD ee (25) 
This equation is multiplied by = Considering (23), the 


following is obtained 


dA _ A dk 
d kak ...(26) 
Wave number conservation (Hutahaean (2024)), 
A 
dk(h+3) pe 
dg 
This equation is expressed as, 
dA 2 ( ‘) dk 3 dh 
dé k dé “dé 
Substituting the left side with (26) 
A dk _ (h N g 
2kd? k 2/dě “dě 
Therefore, 
dk 4k dh 
z Ghee ae ...(27) 


The equation (25) is used to calculate changes in wave 
number k. To compute changes in wave amplitude, 
equation (26) should be utilized. Equation (26) incorporates 
breaking characteristics through A(k, A), where the values 
of k and A determine its behavior. 

As water depth decreases, both k and A increase, causing À 
to decrease continuously until it reaches zero, signifying the 
onset of breaking. Post breaking, A turns negative, resulting 
in a gradual reduction of wave amplitude until it diminishes 
completely. 

For a wave moving from point € to point ë + 6€, the change 
in k and A is 


dk 
Resse = kg + ô ae 
dA 
Azyse = Az + ô aE 
The equation for calculating changes in G is formulated 


using the energy conservation equation which can be 
written as. 
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dG 1dk 


G 2k 

Equation ini Is integrated into,, 
5485 dG 1 (5+5 ak 

Í go L k 

1 

5 (In kesoz — In ke) 


1 
en Gg—-5(In Keyge—In kg) 


In Gë = ]n Gz m= 
Gz+sg = 


8.2. Outcome of the shoaling-breaking model. 

The shoaling-breaking model examines waves with a period 
of 8 seconds and a significant wave height, H = 2.72 m 
(See Table 4), propagating over a bottom slope 


characterized by ~ = —0.02. 


12 


H &L(m) 


0 5 10 
water depth h (m) 


wave height H wave length L 


Fig (2) Shoaling-breaking analysis 
The calculation employs parameters 0 = 1.75 and € = 
0.01. Fig (2) illustrates the outcomes of the shoaling and 
breaking analysis. 


The shoaling-breaking analysis results are depicted in 
Figure (2), where the breaker height Hp, =3.165 m, breaker 


depth h, =6.387 m, and the breaker depth index = 
b 
0.496. This differs notably from McCowan's (1894) 


xk H 
criterion, where oa = 0.78. 
b 


According to Komar and Gaughan (1972), the breaking 
wave height is given by: 


H, = 0.39 g /5(T)H2)/s 
Using Ty = 8.0 sec., Hp = 2.72 m, g = 9.81 m/sec” Hp = 
3.15 m. The breaking wave height from the model closely 


aligns with the Komar and Gaughan (1972) equation. 


Next, we will examine the results of the shoaling-breaking 
model across several wave periods, focusing on the 
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maximum wave height in each period as detailed in Table 
(4). The calculations maintain consistent parameters 0 = 
1.75 and £ = 0.01. The model outcomes are summarized in 
Table (5). 

Table (5) Shoaling breaking outcomes 


T | Ho H, hp H, | Ho-xe 
(sec) | (m) (m) (m) hp (m) 
6 | 153 | 1.781 | 3.593 | 0.496 | 1.772 
7 | 2.082 | 2.426 | 4.896 | 0.496 | 2.411 
8 | 2.72 | 3.165 | 6.387 | 0.496 | 3.15 
9 | 3.444 | 3.995 | 8.06 | 0.496 | 3.988 
10 | 4.252 | 4.93 | 9.95 | 0.496 | 4.923 
11 | 5.146 | 5.945 | 11.996 | 0.496 | 5.958 
12 | 6.124 | 7.081 | 14.289 | 0.496 | 7.091 
13 | 7.186 | 8.331 | 16.811 | 0.496 | 8.321 
14] 8.334] 9.662 | 19.499] 0.496 | 9.65 
15 | 9.568 | 11.081 | 22.363 | 0.496 | 11.079 


Based on the results of the shoaling-breaking model, 2 
b 


: i bi Hp 
remains constant across all wave periods, specifically a 
b 


0.496. The breaker height H, closely matches the Komar- 
Gaughan (1972) breaker height H,_x¢, with a minimal 
difference typically in the second decimal place. 


: H fest 
To achieve = closer to McCowan's (1894) criterion, the 
b 


model is executed using a deep water coefficient 0 = 
1.182. 


Table (6) presents the results of the shoaling-breaking 
analysis 


T Ho Hp hp H, Hy_Ke 
(sec) | (m) (m) (m) hp (m) 
6 1.528 | 1.779 2.28 0.781 1.77 
7 2.08 2.42 3.101 0.78 2.409 
8 2.718 | 3.152 4.038 0.78 3.148 
9 3.44 3.988 5.111 0.78 3.984 
10 | 4.246 | 4.932 6.32 0.78 4.918 
11 5.138 | 5.965 7.645 0.78 5.951 
12 6.116 | 7.085 9.08 0.78 7.083 
| 13 | 7.178 | 8.313 | 10.652 | 0.78 | 8.313 | 
| 14 | 8.324] 9.65 | 12.368 | 0.78 | 9.641 | 
15 9.556 | 11.073 | 14.192 | 0.78 | 11.068 
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It shows a reduction in Hy and H, with a relatively minor 
decrease, while there is a significant decrease in breaker 


depth hp, resulting ae = 0.78, aligning with McCowan's 
b 


criterion. 


The initial criterion for determining @ is based on wave 
dk(h+z) 
dx 


number conservation, where = 0, implying k(h + 


z) = constant. This condition is crucial for solving the 
Laplace equation using separation of variables. For deep 


water, tanhk (h + 2) = tanh kh ~ 1., leading to k (h + 
$) = kh = Om, where tanh 8r ~ 1. Typically, 0 > 1.75 
satisfies this condition. However, from the shoaling- 
breaking model results, k (h + $) or = kh does not need to 


be exactly 1. With @=1.182, we find tanh@z = 
0.988104, which is slightly less than 1 but still close. 


To achieve ab = 0.8, one could consider @ = 1.158. 
b 


Ti . H 
However, determining the most suitable 0 and a values 
b 


requires additional data, including laboratory and analytical 
results. 


By adhering to McCowan's breaker depth index (1894), the 
deep water depth is 


1.1822 Ao 
ho = E 
ho is deep water depth , Ay is deep water wave amplitude 


and Ky is deep water wave number. 


8.2 Refraction-Diffraction Model 

The shoaling-breaking analysis discussed focuses on wave 
transformation along the direction of wave propagation, 
specifically along the € axis. Refraction-diffraction 
analysis, on the other hand, examines wave transformations 
within a plane, typically the (x,y) plane. The changes in 
wave number k and wave amplitude A along the horizontal 
axis x, as depicted in Figure 1, are related to the variables & 


and x. 

dk dk 

ae = ae a 
dA dA 

ax = a a 


With these two equations, the values of k,45, and Ay45, 
can be calculated. By knowing the value of ky+sx the value 
of Gy45,can be calculated. 


1 
G see em Gx-30n Kx45x—In ky) 
x+6x 


Wave direction at point x + ôx, calculated by Equation, 
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u 
a = atan (-) 
dv _ du 
da “ax dx 
dx u2 + v2 


In the following section, we present the results of the 
refraction-diffraction model applied to submerged island 
bathymetry (Fig. 3). The model considers a wave with a 
period T = 8.0 sec and amplitude H = 2.4 m, incident at an 
a =0°relative to the horizontal -x. The calculation 
parameters used are 0 = 1.182 and £ = 0.01. 


Figures 4 and 5 depict the results of the refraction- 
diffraction model. These figures demonstrate that the model 
effectively simulates both refraction-diffraction and 
shoaling-breaking phenomena. 


Incoming wave 
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Fig (3) Submerged island bathymetric contour. 
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Fig (4) Wave height contour on the submerged island 
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gs 


Fig (5) 3-D Wave height, on submerged island 


IX. CONCLUSION 


Increasing the accuracy of KFSBC integration leads to 
improvements in both critical wave steepness and breaking 
wave steepness. Moving from zero-order integration 
accuracy to second-order accuracy results in significant 
changes in both critical wave steepness and breaking wave 
steepness. However, the transition from second-order to 
third-order accuracy shows minimal changes in these 
parameters, indicating convergence of integration. 


The first conclusion drawn is that the integration method 
employed is correct and suitable for the analysis. 
Furthermore, the convergence of integration results 
suggests that third-order accuracy can be considered 
sufficient, with the resulting wave amplitude function 
showing excellent performance. 


The critical wave steepness and breaking wave steepness 
obtained represent final values and can serve as reliable 
references. Similarly, the resulting breaking wave height is 
accurate. 


Although analytical equations for calculating breaker depth 
or breaker depth index were not derived in this research, the 
shoaling model provides estimates of breaker depth. There 
exists a correlation between deep water depth and breaker 
depth: greater deep water depth corresponds to greater 
breaker depth, and vice versa. Certainty about breaker depth 
implies certainty about deep water depth, and vice versa. 
Currently, the breaker depth index derived from empirical 
research by previous investigators serves as an approach for 
estimating deep water depth based on breaker depth. 
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